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The general nuclear contact matrices are defined, taking into consideration all partial waves and 
finite-range interactions, extending Tan’s work for the zero range model. The properties of these 
matrices are discussed and the relations between the contacts and the one-nucleon and two-nucleon 
momentum distributions are derived. Using these relations, a new asymptotic connection between 
the one-nucleon and two-nucleon momentum distributions, describing the two-body short-range cor¬ 
relations in nuclei, is obtained. Using available numerical data, we extract few connections between 
the different contacts and verify their relations to the momentum distributions. The numerical 
data also allows us to identify the main nucleon momentum range affected by two-body short-range 
correlations. Utilizing these relations and the numerical data, we also verify a previous independent 
prediction connecting between the Levinger constant and the contacts. This work provides an im¬ 
portant indication for the relevance of the contact formalism to nuclear systems, and should open 
the path for revealing more useful relations between the contacts and interesting quantities of nuclei 
and nuclear matter. 
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I. INTRODUCTION 

Recently, a new variable, called the contact, was de¬ 
fined by Tan [1, 2] for a two component Fermi gas in¬ 
teracting via short-range forces. The contact measures 
the probability to find two unlike fermions close to each 
other. In a series of theorems, called Tan’s relations, 
many other properties of the system, such as its energy, 
pressure, and momentum distribution, were connected 
to the contact. Tan assumes that the range of the in¬ 
teraction is much smaller than the scattering length and 
the averaged distance between the fermions. Following 
these theoretical predictions, several experiments were 
conducted, verifying Tan’s relations in ultracold atomic 
systems consisting of [3, 4] and ®Li [5-7] atoms. 

Nuclear systems differ from these ultracold atomic sys¬ 
tems in many aspects. First, the nucleons are not two- 
component fermions. Second, while in the atomic sys¬ 
tems the strength of the interaction between the atoms 
and the density can be changed easily, such that Tan’s 
assumptions are satisfied, in nuclear physics it cannot 
be done. In nuclear systems, the s-wave spin-singlet and 
spin-triplet scattering lengths are about —20 fm and 5.38 
fm, respectively, and the average distance between two 
adjacent nucleons is about 2.4 fm. The interaction range 
of the long range part of the nuclear potential, which is 
governed by the pion exchange Yukawa force, is about 
= h/rriT^c ~ 1.4 fm. Thus, in nuclear physics the in¬ 
teraction range is only slightly smaller than the average 
distance between two particles and the scattering length. 
Consequently, some changes are to be done in order to 
generalize Tan’s relations to nuclear systems. 

Considering a two-component Fermi gas that obeys 
Tan’s assumptions, the high momentum tail of the mo¬ 
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mentum distribution is connected to the contact through 
the relation, n„{k) —^ C/k‘^ as A: —)■ oo, where n„{k) is 
the momentum distribution of fermions with spin a, and 
C is the contact. In nuclear physics, the high-momentum 
part of the nucleon’s momentum distribution is one of the 
main tools for studying short range correlations (SRCs) 
between nucleons. The main focus in current studies of 
two-body SRCs (see e.g. [8-12]) is around the momen¬ 
tum range 1.5 fm“^ < fc < 3 fm“^. In few of these stud¬ 
ies it is claimed that higher momentum is affected also 
by 3-body correlations [13]. In this momentum range, a 
dominance of neutron-proton {np) correlated pairs was 
observed in electron scattering experiments [10]. This 
np dominance is usually explained by the contribution of 
the tensor force, which affects only spin-triplet np pairs. 
Another observation is that the correlated pairs usually 
have high relative momentum and low center of mass 
momentum, i.e. they move approximatelly back-to-back. 
Generalizing Tan’s relation between the high momentum 
tail and the contact to nuclear systems, should help in 
understanding more properties of SRCs in nuclei. 

In a previous paper [14], we have suggested that it 
might be fruitful to use the contact formalism in nuclear 
systems. There we have defined the neutron-proton s- 
wave nuclear contacts and evaluated their average value 
by relating them to the Levinger constant of the photoab¬ 
sorption process. In this work we generalize the definition 
of the nuclear contacts from s-wave to all partial waves. 
We also consider finite-range interactions instead of zero- 
range. The result is the matrices of nuclear contacts. We 
discuss the properties of these matrices, and use our gen¬ 
eralized contact formalism to relate the nuclear contacts 
to the one-nucleon and two-nucleon momentum distribu¬ 
tions in nuclei. Doing so, we find an asymptotic relation 
between these two distributions which is relevant to the 
study of SRCs in nuclei. This relation is verified by avail¬ 
able numerical data. Further analysis of the numerical 
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data and its implications to the contact formalism are 
also presented. 

In this paper we focus on the two-body contacts and 
on two-body correlations, postponing the discussion on 
three-body effects to future publications. 

II. THE MATRICES OF NUCLEAR CONTACTS 

Consider a two-component Fermi gas that obeys Tan’s 
assumptions. In such a gas, when a spin-up particle i 
gets close to a spin-down particle j, the many-body wave 
function can be factorized into a product of an asymp¬ 
totic pair wave function ipirij), rij = u — r^, and a 
function A, also called the regular part of describing 
the residual A — 2 particle system and the pair’s center 
of mass Rij = {vi + rj)/2 motion [1, 15], 

^ ■ ( 1 ) 

Due to the suppression of higher partial waves in these 
systems, the asymptotic pair wave function will be pre¬ 
dominantly an s-wave. In particular, in the zero-range 
model [16] the pair wave function is given by v? = 
(l/r^ — 1/a), where a is the scattering length. 

The contact C is then defined by [1, 15] 

C' = 16^2^tf(^l^)> (2) 

where 



X A^ (^Rij^ ■ A (J?j) 

and Njij is the number of possible spin up - spin down 
pairs. 

In nuclear physics, we have four-component fermions, 
which are the protons and neutrons with their spin being 
either up or down. Moreover, the assumption of a zero- 
range s-wave interaction is not accurate for nuclei. As a 
result, few changes must be made in order to generalize 
the contact formalism to study nuclear systems. 

A nucleus can be described by a wave function with 
total angular momentum J and projection M. We will 
assume that when particle i gets close to particle j, the 
wave function is still factorized but the pair wave func¬ 
tion depends on the total spin of the pair S 2 , and its 
angular momentum quantum number £2 (with respect to 
the relative coordinate r^) which are coupled to create 
the total pair angular momentum j 2 and projection m 2 . 
The asymptotic form of the wave function is then given 

by 

^ ^ ■ ( 4 ) 

Tij —>-U 

a 

Here the index ij corresponds to one of the three par¬ 
ticle pairs: proton-proton (pp), neutron-neutron (nn) or 


neutron-proton {np). We note that due to symmetry the 
asymptotic functions are invariant under same particle 
permutations. The sum over a denotes a sum over the 
four quantum numbers {s 2 ,£ 2 , j 2 ,'m 2 ). 

= E (j2m2J^-2M^_2lJM)Alf 

J A-2,Ma-2 

( 5 ) 

Here, Ja -2 and Ma -2 are the angular momentum quan¬ 
tum numbers with respect to Ja- 2 +L 2 ,cm, where Ja -2 
is the total angular momentum of the residual (A — 2) 
particles and L 2 ,cm is the spatial angular momentum 

with respect to Rij. is a set of func¬ 

tions with angular momentum quantum numbers Ja -2 
and Ma- 2 , which depends also on the numbers S 2 ,l' 2 )j 2 - 

^ , ( 6 ) 

where Xs 2 )j-a ^be two-body spin function, and 

For clarity, 

when angular momentum indices are written without 
any brackets they denote the relevant angular momen¬ 
tum quantum numbers of the function. When the indices 
are in curly brackets, it means that the function depends 
on this numbers but they do not denote the angular mo¬ 
mentum of the function. When two indices are inside 
round brackets, it means that the angular momentum of 
the function is created by a coupling of these two indices. 

The only assumption we make regarding the set of 
functions Wfj} is that they do not depend on the specific 
nuclei or its total angular momentum J and M. This is 
a reasonable assumption, because when two particles are 
very close they interact with each other regardless to the 
background of the A — 2 particle system. Doing so, we 
no longer use the s-wave or the zero-range assumptions. 

Since the A“ functions are not generally orthogonal 
for different a, we are led to define matrices of nuclear 
contacts in the following way: 

C'“^(JM) = 16^2^y(A“lAf^.). (7) 

As before, ij stands for one of the pairs: pp, nn or np, 
Nij is the number of ij pairs, and a and /3 are the matrix 
indices. We also denote a = {sa^£aTja,'>na) and j5 = 
{sp,£i 3 ,ji 3 ,mp). One can see that if nia m^, then 
C^j^iyJM) = 0, but it is not generally true for j 2 , S 2 or 

£ 2 - For spherical nuclei (J = 0) we do get Cff{JM) = 
0 if ja 7^ jp- For pp and nn pairs, Pauli’s exclusion 
principle tells us that unless Sq -I- £a is even, we have 
= 0, so = 0 if or sp + £p 

are odd. Moreover, if 'I' is the ground state of the nucleus, 
or any eigenstate of the nuclear Hamiltonian, then 'k has 
a defined parity. has a parity of (—1)^“, so it dictates 

the parity of A“-. Thus, C'“^(JM) = 0 for a and /? such 
that £a and £p have different parities. 

Since the projection M is usually unknown in experi¬ 
ments, it is useful to define the averaged nuclear contacts: 

M 
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According to this definition, we have three matrices of 
averaged contacts, one for each kind of nucleon-nucleon 
pair. We note that the averaged contacts still depend on 
J, but we will not write it explictly. Using Clebsch Gor- 
dan identities one can prove that if rria ^ mp or ja ^ jp, 
then = 0, and also that the averaged contacts are in¬ 
dependent of rria and mp. The averaged contacts inherit 
the properties of the non-averaged contacts C^^^{JM) 
regarding parity and Pauli’s principle. 

Concluding, for a given a, the relevant /3’s such that 
can be different from zero must obey jp = ja and 
mp = ma- Since, S 2 = 0,1 there are four ( 52 ,^ 2 ) pairs 
that can create a given ^ 0: {Q,ja), (1, ja), (1, ja - 1) 
and (1, ja + l). The first two options have the same parity 
of £2 and the last two have the opposite parity. For j„ = 0 
we have only two possible (s 2 ,£ 2 ) pairs: (0,0) and (1,1), 
which have different parity of £ 2 - Thus, in general the 
matrices are built from 2x2 blocks, except for the 
two 1x1 blocks associated with the j 2 = 0 case. Each 
block has a well defined j 2 ,m 2 values. For any j 2 0 
there are two blocks, one with ( 52 ,^ 2 ) = (0, j 2 ), (1, j 2 ) 
and the other with (s 2 ,l' 2 ) = (1, j 2 —1), (1, j 2 + l)- For pp 
and nn pairs, Pauli’s principle dictates that any matrix 
element with an odd S 2 + £2 is zero, so some of the 2x2 
blocks are reduced into two 1x1 blocks. 

In a previous paper [14] we have defined the s-wave 
nuclear contacts, C-^{JM), for S 2 = 0,1. The definition 
there was slightly different from the current one, as the 
two-body spin functions were included into the regular 
(A—2) particle function A“ . In our current definition, we 
have four diagonal s-wave contacts and 

where aoo = (s 2 = 0,^2 = 0 ,j 2 = 0 ,TO 2 = 0), ai^ = 
(s 2 = 1, £2 = 0,j2 = 1,TO2 = p), and p = —1,0,1. The 
relations between the two definitions are 

Clp°{JM) = (9) 


1 

Clp^{JM)= Y, ( 10 ) 

ii=-i 

Averaging over M and using the fact that the averaged 
contacts are independent of m 2 we get 

^S2—0 _ ^aooQ^oo (11) 


H (12) 

We also note that the previously defined s-wave contacts, 
C'-j{JM), are actually independent of M. Thus, also 
C^ooo.oo{JM) and Ei=-i C'“^“''^(JM) are independent 
of M. 


III. MOMENTUM DISTRIBUTIONS 
A. The two-nucleon momentum distribution 


In the following we will utilize the above generalized 
contact formalism to find a relation between the two- 
nucleon momentum distribution and the nuclear con¬ 
tacts. 

Let’s denote by f£^\k + K/2,—k + K/2) the density 
probability to find a pair of nucleons, ij G {pp, nn,pn}, 
with any particle of type i with momentum k + K/2 and 
any particle of type j with momentum —k + iT/2. J 
and M are the angular momentum quantum numbers of 
the nuclear wave function 'k. Working in the momentum 
space 

^(fci,...,fc^)= / (13) 


and we can write 


PP[k + KI2,-k + KI2) = N,, 


n 

m^i,j 


d kYn 

{2ttY 


^(fci, ...,ki = k + K/2, ...,kj = —k + K/2, ...,kA) 


(14) 


where A is the number of nucleons, N^j is the number of 
ij pairs, and we notice that is normalized in such 

away that / ^ ^ = Wj ■ 

In the limit A: —>■ 00 the main contribution to 
comes from the asymptotic rij —> 0 part of the wave 
function, given in Eq. (4). All other terms will cancel 
each other due to the fast oscillating exp(ifc • r^) factor. 
Substituting ^ into Eq. (14), and using Eq. (4) we get 



X exp(ife • -I- iK ■ Rij + ^ ikn ■ rn)\^. (15) 


We will define now 

^)f'(fc)= J ^f/,^{k + K/2,-k + K/2). (16) 

is the density probability to find an ij pair with 
relative momentum fc, and it obeys the normalization 
condition J Fd^(k) = Nij. We can now substitute 
the asymptotic form of Eq. (15), into the defini¬ 
tion of Pd^. In the resulting expression we can separate 
the integration over from the rest of the coordinates. 
Using the notation 

£p%{k) = J d^rpf^ir) exp{ik ■ r) (17) 
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and 


= y n d^'r’nd^RijA'^jexp{iK-Rij+ ^ iK-Vn) , 

(18) 

we get 




a,(3 


m^i,j 


(Pkm d^K 


(27r)3 (27r)3 


( 19 ) 


where iiTps is the center of mass momentum of the ps pair. 
Substituting into np^{k), we see that since Ap^ is regu¬ 
lar, Ap^ will be significant only if |fc -|- fe^l <C fc. It means 
that kg « —k so k — kg m 2k. Substituting into Eq. 
(23), we get summations over s, s' ^ p. The contribution 
of the s, s' element, for s ^ s', will be significant only 
for kg « kg' « —k. In this case k, kg, kg' —>■ oo together, 
which is clearly a three body effect and we expect it to 
be less important [17]. Therefore we are left only with 
the diagonal elements and get 


Noting the equality 


n d kjyi d K raf 

(27r)3 (27r)3 b b 

m^i,j ^ ^ ' 



( 20 ) 

we obtain the following asymptotic k ^ oo expression 
for the two nucleon momentum distribution. 


F/,^{k)=Y,ptJ{k)p^^{k) 


IGtt^ 


( 21 ) 


Here we have used the definition of the contacts from Eq. 
(7). Averaging over M, we get the asymptotic relation 

F,,(fe) = ^<p“l(fc)^f^.(fc)^, (22) 

a,j3 




s^p Q!,/3 


(27r)3 (27r)3 


X AfgiKpg, {k,},^p,g)AUKpg, {k,},^p^g). 


(25) 


We will now divide the sum i^^to a sum over pro¬ 

tons and a sum over neutrons X)p'^p + X)n- Since the 
asymptotic functions App, and iPpp, are the same for all 
pp' pairs we can take them out of the sum. The same 
holds for the np pairs. As a result we get 


<^(fc) =J2pfp{k)p^pp{k)Z(Z-l){A;p\APp) 

Oi,j3 


xJ2p;l{k)p^p„ik)NZ{A^JA^pJ. (26) 

Ct,0 


where = (2 J -|- 1) , and are the av¬ 

eraged contacts defined in Eq. (8). Like C“^, also 
depends implicitly on J. 


B. The one-nucleon momentum distribution 


We would like now to connect the nuclear contacts 
also to the one-nucleon momentum distributions. The 
following derivation is based on Tan’s derivation for the 
two-body case in atomic systems [1]. We will first fo¬ 
cus on the proton’s momentum distribution rip^^k). 
Normalized to the number of protons in the system Z, 

/ (Cf'^p^^(^) = "■p^ is i’y 



^{ki,...,kp = k,...,kA) , 

(23) 


where p is any proton. 

In the k —> oo limit the main contribution to rip^ 
emerges from the asymptotic parts of the wave function, 
i.e. from rpg = jr^ — r^j —>■ 0, for any particle s A 
being proton or neutron. In this limit 


^(fci, ...,kp = k, ..., kA) = J2J2 ((^ - fc«)/2) 

S^p Oi 

^ ^ps i^ps = k kg, {kj}j^p^g ), 


Here N is the number of neutrons in the system. Using 
the definition of the contacts, Eq. (7), we see that for 
k ^ oo 

rip {k) = Ppl{k)plp{k} 

+ (27) 

a,/3 

Averaging over M we further obtain the relation between 
the averaged contacts and the averaged protons’ momen¬ 
tum distribution np{k) = (2J -|- 1)“^ '^p^(^) 

k —>■ C!o: 

2(7“^ 

ripif^) = '^Ppli^)‘Pppi^)j^ 

+ J2^pnik)p^pnik)Y^. (28) 

a, 13 

We note that np{k) still depends on J. Similarly, for the 
neutrons: 


n„(fc)=^^“t(fc)<pL(fe)^ 

+ Y.^pn{k)pluik)^. 

Ol,0 


(24) 


( 29 ) 
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Comparing Eqs. (28) and (29) to Eq. (22), we can 
see that for k —> oo there is a simple relation between 
the one-nucleon and the two-nucleon momentum distri¬ 
butions: 

np{k) = 2Fpp(k) + Fpn{k) (30) 

n-aik) = 2Fnn{k) + Fpn{k). (31) 

These connections seem intuitive if we assume that a nu¬ 
cleon will have high momentum k only if there is another 
nucleon close to it with opposite momentum —k. In this 
case, if we find a proton with high momentum k we know 
that we will find close to it a neutron or a proton, that 
is a correlated pp or np pair with relative momentum k. 
Notice that the factor of 2 before Fpp and in Eqs. 
(30), (31), can be also explained in this picture by the 
fact that for example a pp pair with momenta (—fc, k) 
has a relative momentum —k even though there is a pro¬ 
ton with momentum k in this pair. It means that such 
a pair will be counted for np(k) but not for Fpp(k) and 
the factor of 2 takes it into consideration. 

These relations emphasize the importance of the two- 
body correlations to the high momentum one-nucleon 
distribution. As mentioned before, the picture of short- 
range correlated pairs of nucleons with back-to-back mo¬ 
mentum is one of the main features of SRCs in nuclei, 
and the above relations between the one-nucleon and 
two-nucleon momentum distributions give a theoretical 
support to this picture. 

We also note here that similar derivations can be done 
easily for atomic systems consisting of two-component 
fermions, denoted by f and The one-body high 
momentum distribution is already known and given by 
n^(fc) = n^(fe) = Cjk'^. Adjusting the above derivation 
for the two-nucleon momentum distribution to atomic 
systems will produce an identical relation between the 
two-body momentum distribution, E^^(fc), describing the 
probability to find an fj, pair with high relative momen¬ 
tum, and the atomic contact. Explicitly, Ef^(fc) = C/k'^. 
As a result we find that n^(fe) = ni{k) = E|^(fc) for 
high momentum k. This relation tells us that also in 
the ultracold atomic systems the correlated pairs have 
back-to-back momentum, like in nuclear systems. 


IV. ANALYSIS OF NUMERICAL DATA 

A. Momentum distributions 

In order to check the validity of our results in actual 
nuclear systems, we turn now to compare our theoretical 
predictions to available numerical data. To this end, we 
will use numerical data of one-nucleon and two-nucleon 
momentum distributions calculated by Wiringa et al. [8] , 
using the Variational Monte Carlo method (VMC), for 
nuclei with A < 12. In this VMC results, the calculation 



FIG. 1. (Color online) The ratio {2Fpp+Fp„)/np for different 
nuclei. The numerical data is taken from Ref. [8]. Red line 
- ‘*He, green line - ®He, cyan line - ®He, black line - ®Li, blue 
line - ®Be, and pink line - The dashed red line is the 
reference y = 1. 



FIG. 2. (Color online) The ratio (2F’„„ + Fpn)ln„ for the 
non-symmetric nuclei in the numerical data of Ref. [8]. Blue 
line - ®He, and green line - ®He. The dashed red line is the 
reference J/ = 1. 


of both one-nucleon and two-nucleon momentum distri¬ 
butions were done for nuclei in their ground state. Con¬ 
sequently, the following analysis is limited to the nuclear 
ground state. 

First we check the relation between the one-nucleon 
and the two-nucleon momentum distributions, Eqs. (30), 
(31). In Fig. 1 the ratio between 2Fpp + Fpn and Up is 
presented for various nuclei. We can see that for k —> oo 
the two quantities coincide and our prediction (30) is 
indeed satisfied. In Fig. 2 we present the ratio between 
‘2Fnn + Fpn and We show only the results for non- 
symmetric nucei, because for symmetric nuclei there is no 
difference between protons and neutrons in the numerical 
VMC data. We can see that also here, the ratio (2F„„ + 
Fpn)/nn —> 1 as fc —> oo and our prediction (31) is 
satisfied. This result is obtained for all available nuclei: 
"^He, ®He, ®He, ®Li, ®Be, and for both protons and 
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neutrons. For all these nuclei the momentum relations 
hold for 4 fm“^ < k < 5 fm“^. 

The correspondence between our predictions, derived 
using the contact formalism, and the numerical data is a 
good indication for the relevance of the contact formalism 
to nuclear systems. We also learn here that the approx¬ 
imations made in the above theoretical derivations for 
k —>■ oo are valid for 4 fm“^ < k < 5 fm“^. This is the 
first indication for the momentum range which is relevant 
to the contact formalism in nuclear systems. Moreover, 
as we mentioned before, in current studies of SRCs in 
nuclei this momentum range of fc > 4 fm“^ is believed 
to be affected by three-body correlations. As explained, 
Eqs. (30) and (31) are suppose to be satisfied when the 
two-body correlations are the only significant correlations 
and every high momentum nucleon has a sinlge nucleon 
near it with back-to-back momentum. It means that ac¬ 
cording to this numerical data the momentum range of 
4 fm“^ < fc < 5 fm“^ is affected almost exclusively by 
two-body SRCs while three-body SRCs are negligible, 
and that in this momentum range the picture of back-to- 
back short-range correlated pairs is accurate. 

We note that this momentum range of 4 fm“^ < fc < 
5fm“^ might be model dependent, and it should be veri¬ 
fied using other numerical methods, and different nuclear 
potentials. It should also be mentioned that the VMC 
method utilize two and three-body Jastrow correlations 
in the nuclear wave function. 

Hen et al. [18] also discuss the possibility that the 
contact formalism is relevant in nuclear physics. In their 
work, they present an experimental measurement of a 
k~^ behavior in the proton momentum distribution in 
the deuteron for 1.6 fm“^ < fc < 3.2 fm“^. They also 
claim that the k~^ behavior exists in heavier nuclei in 
the same momentum range. As mentioned before, one 
of the results of the contact formalism in atomic systems 
is the k~^ tail in the momentum distribution, but this 
behavior is a direct consequence of the zero-range model. 
In nuclear systems this model is not accurate, so we can 
only expect a high momentum tail universal to all nuclei, 
but not a k~^ behavior. We also note that in the nu¬ 
merical VMC data there is no k~'^ tail for nuclei heavier 
than the deuteron. Moreover, we have found here that 
the relevant momentum range for the contact formalism 
in nuclei is 4 fm“ ^ < fc < 5 fm“ ^, which is higher than 
the momentum range discussed by Hen et al. 


B. The pp and nn contacts along the nuclear chart 

We continue now by examining the ratio between 
Fpp{k) and Fnn{k) in the same nuclei. In the VMC re¬ 
sults, this ratio equals 1 for all fc for symmetric nuclei 
{N = Z). Therefore, we are left with the available non- 
symmetric nuclei ®He and ®He. The relevant results are 
shown in Fig. 3. 

We can see that for 4 fm“^ < fc < 5 fm“^ the ra¬ 
tio is approximately constant. Inspecting Eq. (22), 



FIG. 3. (Color online) The ratio between Fpp and F„„ in 
the same nuclei for the available non-symmetric nuclei in [8]. 
Blue line - ®He, and red line - ®He. The dashed blue and red 
lines indicate the value of Z/N in ®He, and ®He, respectively. 

we see that the only way for this ratio to be constant 
is that (i) only pairs in a,/3 states with the same fc- 
dependence of contribute significantly to both Fpp 

and Fnn, and (ii) both pp and nn pairs have the same 
fc-dependence. It is reasonable to assume that the s- 
wave contacts are the most significant contacts. For pp 
and nn pairs the only possible non-zero s-wave contact 
is C"™"™ where aoo = (s 2 = 0,^2 = 0 ,j 2 = 0 ,m 2 = 0). 
This point can be verified numerically through analysis of 
the angular dependence of the momentum distributions. 
If the s-wave contact is indeed dominant we expect to 
see no angular-dependence. If we further assume that 

which seems reasonable from 
isospin symmetry, then the ratio between Fpp and Fnn 
for large momentum equals to the ratio between Cp°°°‘°° 
and (7“™“™. 

We can also see in Fig. 3, that for the two relevant 
nuclei the ratio Fpp/Fnn is close to the ratio Z/N between 
the number of protons and neutrons in the nucleus. If this 
relation turns out to be true in general along the nuclear 
chart, it means that for a nucleus X in its ground state, 
the most significant pp and nn contacts are and 

and their ratio is given by 


C'aoo«oo(X) 

z{x) 

(32) 


N{X)' 


aooQ^oo 
rnn ' 

(33) 


Here Z{X) {N{X)) is the number of protons (neutrons) 
in the nucleus X. This result is surprising because one 
might think that the ratio (32) should scale as the ratio 
between the total number of pp pairs and the number of 
nn pairs in the nucleus, i.e. Z^/N"^. The above result 
tells us that the number of correlated pp and nn pairs in 
nuclei goes like Z and N, respectively. If we check the 
ratio between Fpp or Fnn and Fpn, no plateau is observed. 
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FIG. 4. (Color online) The ratio between Fpp{X) and 
F'ppC^He) for the available nuclei X in the numerical data of 
Ref. [8]. Blue line - ®He, red line - ®He, green line - ®Li, black 
line - ®Be, and pink line - ^'^B. 


We can also examine the ratio between Fpp of nucleus 
X and Fpp of another nucleus Y. The results are pre¬ 
sented in Fig. 4, where all the available nuclei are com¬ 
pared to ^He. 

Here again we see flattening for 4 fm“^ < fc < 5 fm“^. 
This behavior supports the claim that only one contact 
contributes significantly to Fpp^ and so the value of this 
ratio is just the value of the ratio of this pp contact in 
the two different nuclei (see Eq. (22)). The constant be¬ 
havior also supports the assumption that the pair wave 
functions Pp^ are universal along the nuclear chart, be¬ 
cause that way the fc-dependence indeed vanishes. The 
average values of this ratio for 4 fm“^ < k < 5 fm“^ 
are presented in table I and compared to the ratio be¬ 
tween the number of protons in the relevant nuclei and 
the number of protons in "^He. We can see that the two 
ratios are approximately equal for the different nuclei. 
If the most significant pp contact is the s-wave contact 
then we we can deduce that for nuclei X and Y 
in their ground state: 

C^ooo^oo(Y) Z(Y)' ^ ’ 


For Fnn similar results are observed, therefore we can 
also deduce that 


Caooaoo(X) N{X) 

C'aooaoo(y) ~ ' 


(35) 


These relations support the claim that the number of 
correlated pp and nn pairs in nuclei is proportional to Z 
and TV, respectively. 


C. The pn contacts and the Levinger constant 

So far we have studied the properties of the pp and 
nn contacts, now we turn to study the pn contacts. The 


X (Fpp(X)/Fpp(‘^He)) Z(X)/Z(^He) 


'^He 

0.94±0.01(lcr) 

1 

®He 

0.90±0.01(lcr) 

1 

6Li 

1.09±0.02(lcr) 

1.5 

®Be 

2.31 ± 0.07(lcr) 

2 

10b 

2.91 ± 0.06(lcr) 

2.5 


TABLE I. The averaged value of the ratio between Fpp{X) 
and Fpp ('^He) for 4 fm“^ k < 5 fm“^ for all the available 
nuclei in the numerical data of Ref. [8]. 



FIG. 5. (Golor online) The ratio between Fp„{X) and np(^H) 
for the available nuclei X in the numerical data of Ref. [8]. 
Blue line - '^He, red line - ®He, green line - ®He, black line - 
®Li, pink line - ®Be, and cyan line - ^°B. 


pn contacts might be the most interesting ones because 
of the dominance of correlated pn pairs in nuclear SRCs 
[10]. In order to study the properties of the pn con¬ 
tacts we examine the variation in Fpn between different 
nuclei. As in the pp and nn cases, also in this case we 
shall assume that the s-wave is the most dominant partial 
wave. For a pn pair in an s-wave there are two possible 
spin configurations, spin-singlet and spin-triplet. For the 
deuteron ^H, only the spin triplet is relevant as it is a 
J = 1 state. In Fig. 5 we present the ratio between 
Fpn{X) and np(^H) for the available nuclei in the VMC 
results. 

Once again, a constant behavior is seen for 4 fm“^ < 
fc < 5 fm“^. As mentioned before, we have three equal 
spin triplet s-wave np contacts, and one spin- 

singlet s-wave np contact 0^°°°“°°. Moreover, 
is independent of p. Consequently, we would expect 
to see a plateau in the ratio Fpp(X)/np(^H), if either 
the asymptotic pair wave functions obey the relation 
|<^pn'"P = \Ppn\^ or alternatively if the spin-triplet s- 
wave contacts are dominant. In the first case we can de¬ 
duce from the relations between the contacts and the one- 
nucleon and two-nucleon momentum distributions that 






















asymptotically 


Fp„(X) + 

np(2H) ^ 3C'“^““^“(2H) 

ctr°{x) + ctr^{x) ^ ^ 

c^r'(2H) 

where here we have also used the notation of Eqs. (11) 
and (12). In the second case we get 

Fp^jx) _ ^ c;r^{x) 

np(2 ii-)c“^““i“(2 H) c;r\^Hy ^ ^ 

In a previous paper [14], we have predicted that the ra¬ 
tio between the sum of the two s-wave np contacts of a 
nucleus X in his ground state and the deuteron’s s-wave 
np contact is given by 


c;r°ix) + c;rHx) _^nz 

di ’ 


(38) 


where L is Levinger’s constant that relates, at the high 
energy hand, the photoabsorption cross section of a nu¬ 
cleus to the photoabsorption cross section of the deuteron 
[19]. Analysis of the experimental results [20] suggest 
that the L is approximately a constant along the nuclear 
chart L « 5.50 ± 0.21 [14]. In [14], we have assumed 
that the two s-wave states have the same asymptotic 
pair wave function in small distances, which corresponds 
to the first case above. If we were to assume that only 
the spin-triplet np s-wave is significant, then our result 
would have been 


c;rHx) _^Nz 

A ■ 


(39) 


In any of the two cases, we get the relation 


FpjX) _ NZ 

np(2H) ^ A 


(40) 


that should hold in the high momentum range. For 
this range of high momentum the ratio between Fpn 
and np(^H) is the number of quasideuteron (qd) pairs 
with high relative momentum in the nucleus. In ta¬ 
ble II we present the averaged value of this ratio for 
4 fm“^ < k < 'b fm“^ and its multiplication by A/NZ 
for each nuclei X, which should be equal to L according 
to the above prediction. One can see that the values of 
the multiplied ratio are close to the above value of L for 
all the nuclei and their average value is 5.7 ± 0.7(lcr). 
This value is in a very good agreement with the above 
mentioned value of L. 

Evaluation of Levinger’s constant from the number of 
qd pairs was done by Benhar et al. [21]. In their work, 
they calculate numerically the number of qd pairs in the 
nucleus and extract Levinger’s constant. In our eval¬ 
uation we consider only the qd pairs with high relative 
momentum, which corresponds to small relative distance. 


A 

{Fp^{X)/np(^H)) A/NZ{Fp„iX)/npCH)) 


6.10±0.06(la) 

6.10 ±0.06 

6He 

6.5±0.1(lcr) 

4.88 ±0.08 

®He 

7.82±0.03(la) 

5.21 ± 0.02 

®Li 

7.63±0.04(la) 

5.09 ±0.03 

«Be 

13.25 ±0.08(1(7) 

6.63 ±0.04 

10b 

15.3±0.3(1(t) 

6.1 ±0.1 


TABLE II. The averaged value of the ratio between Fpn{X) 
and np(2H) for 4 fm“^ < fc < 5 fm“^ and its multiplication 
by AjXZ for all the available nuclei in the numerical VMC 
results of [8]. 


Only such qd pairs can be emitted in the photoabsorption 
process, and therefore only they should be considered. 

We have compared here two independent relations be¬ 
tween the np contacts and different properties of nuclei 
(momentum distribution and photoabsorption cross sec¬ 
tion) and obtained a good agreement between the two. 
Doing so, we have also obtained here an established es¬ 
timation for the leading s-wave np contact (s) along the 
nuclear chart for nuclei in their ground state (in units of 
the deuteron’s s-wave np contact). 


V. SUMMARY 


Summing up, we have generalized the contact formal¬ 
ism to nuclear systems and defined a matrix of contacts 
for each particle pair: pp, nn and pn. With this gener¬ 
alization we have taken into consideration both different 
partial waves, and finite-range interaction. We have dis¬ 
cussed the simple properties of the nuclear contacts and 
demonstrated the use of the generalized formalism by re¬ 
lating the contacts to the one-nucleon and two-nucleon 
momentum distributions. As a result we have obtained 
a relation between these two momentum distributions, 
which emphasizes the significant contribution of SRCs 
to the high one-nucleon momentum tail. Using avilable 
VMC numerical data [8], we have verified the above re¬ 
lation and found further relations between the different 
nuclear contacts. Using few of these new relations and 
a previous prediction connecting the pn contacts to the 
Levinger constant, we have calculated Levinger’s con¬ 
stant for the available nuclei and got a good agreement 
with its experimental value. This is an important in¬ 
dication for the relevance of the contact formalism to 
nuclear systems, and might open the path to revealing 
many more interesting relations. We have also learned 
from the numerical data that the relevant momentum 
range for the contact’s approximations in nuclear sys¬ 
tems is 4fm“^ < k < 5fm“^. However, we note that this 
result might be model dependent. The fact that the rela¬ 
tions between the one-nucleon and two-nucleon momen¬ 
tum distribution were satisfied in this momentum range 
teaches us that for such momenta the two-body SRCs, 
rather than three-body SRCs, are dominant. Additional 


















9 


numerical or experimental data for both one-nucleon and 
two-nucleon momentum distributions in more nuclei, also 
in excited states, including angular-dependence is needed 
in order to improve our understanding regarding the 
properties of the nuclear contacts. 
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